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Abstract. 

A short account of the physics of electrical discharges in gases is given from the viewpoint of its 
historical evolution and application to planetary atmospheres. As such it serves as an introduction 
to the papers on particular aspects of electric discharges contained in this issue, in particular in the 
chapters on lightning and the discharges which in the last two decades have been observed to take 
place in Earth’s upper atmosphere. In addition to briefly reviewing the early history of gas discharge 
physics we discuss the main parameters affecting atmospheric discharges like collision frequency, 
mean free path and critical electric field strength. Any discharge current in the atmosphere is clearly 
carried only by electrons. Above the lower boundary of the mesosphere the electrons must be 
considered magnetized with the conductivity becoming a tensor. Moreover, the collisional mean 
free path in the upper atmosphere becomes relatively large which lowers the critical electric field 
there and more easily enables discharges than at lower altitudes. Finally we briefly mention the 
importance of such discharges as sources for wave emission. 


Keywords: atmosphere, quasi-stationary electric fields, discharges, avalanches, sprites-jets, elec- 
trostatic waves 


1. A brief historical overview 


The physics of electrical discharges in gases has a long history (Bowers 1991, 1998) 
reaching back into the 17th century; in 1672 Otto von Guericke reported the production 
of static electric sparks when rubbing sulphur balls with the hand. Three years later the 
astronomer Jean Picard reported that shaking a mercury barometer caused the tube to 
glow. At the end of that century (around 1698) sparks had been extracted from rubbed 
amber and formed the subject of intense investigation by the Comte Charles Francois de 
Cisternay du Fay. Based on Picard’s observation Johann Heinrich Winckler in Leipzig 
invented what we can call a ‘fluorescent tube’ in 1745, and in 1802 Vasilii Petrov in St 
Petersburgh produced electric arcs between two carbon electrodes independently of, but 
almost at the same time as, Sir Humphry Davy (Knight 1992) in London. Davy invented 
the carbon arc lamp, presenting it in 1807 to the Royal Institution, and established the 
subject of arc physics (Davy 1807). For their experiments they both used the chemical 
battery invented in 1799 by Count Alessandro Volta. 

The famous experiments of Benjamin Franklin in 1750 on sparks led him to coin 
the notions of positive and negative charging which survive until today and to invent 
the lightning rod. Mikhail V. Lomonossov in 1743 (Lomonossov 1748) had already sug- 
gested that lightning, and the polar lights (aurorae) as well, would be atmospheric electric 
discharge phenomena. With respect to the similarity between lightning and sparks they 
were both right, while they could not know that auroral discharges are phenomena quite 
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different from spark discharges. Charles-Augustin de Coulomb (1785) discovered dark 
discharges (Gillmor 1971). In the 19th century interest in gas discharges increased due to 
Michael Faraday’s and James Clerk Maxwell’s formulation of the physics of electromag- 
netic phenomena. Faraday himself even experimented with sparks and glow (or coronal) 
discharges, thereby uncovering the dependence of the latter on the pressure in the tube 
(Faraday 1833, 1834). It is interesting to note that such investigations were made possible 
by improved glass blowing techniques the evacuation of glass tubes and the separation of 
gases. 

Experimenting on dilute gases, Johann Hittorff (Hittorff 1879) determined the elec- 
trical conductivity of air and different gases which, in 1900, culminated in Paul Drude’s 
seminal formula for the electrical conductivity of a gas, introducing the concept of colli- 
sion frequency v,. The collision frequency received its physical explanation when Ernest 
Rutherford found the now famous Rutherford collisional cross section, o,, in terms of 
which the collision frequency can be written as v, = no,v, where n is the gas density, 
and v the velocity of the moving particle. Equivalently this allowed the introduction of the 
collisional free flight distance \g = v/v, of a particle between two collisions. When v 
was taken to be the thermal velocity 74, this became the mean free path Amfp = Vth [Ve = 
1/noc. 

Hittorff somewhat later published the first investigations on cathode rays which were 
crucial for the later formulation of atomic physics. Indeed, any deeper understanding of 
gas discharges had to wait until the advent of atomic physics, which became possible only 
via reference to the investigation of gas discharges. William Crookes (Crookes 1879) 
found that cathode rays were deflected by magnetic fields and called them the fourth 
state of matter, claiming that their investigation would lead to the deepest insight into 
the nature of matter. This prediction was in fact true. From the deflection of cathode 
rays in a magnetic field and based on the electrodynamics of charged particle motions in 
electromagnetic fields, Sir Joseph Thomson (Thomson 1897) determined the elementary 
charge-to-mass ratio which led to the discovery of the electron. One year earlier Wilhelm 
Conrad Roentgen (Roentgen 1895) experimenting with cathode rays had discovered X- 
rays. 

The last decade of the 19th century and the first decade of the next were devoted to 
investigations of spectral phenomena with the help of dilute gas tubes, leading on the 
one hand to the determination of the Johann Balmer (Balmer 1885), Walther Ritz (Ritz 
1908), Theodore Lyman (Lyman 1914) and Friedrich Paschen (Paschen 1908) series of 
spectral line, and on the other hand to the precise mapping of the temperature dependence 
of black body radiation, and the energetic of collisional excitation of electrons in atoms 
at specific energies by James Franck and Gustav Hertz (Franck and Hertz 1914). These 
investigations culminated in the formulation of quantum and atomic theory, going hand in 
hand with the development of gas discharge physics. It is thus no surprise that most of the 
Nobel laureates of the first three decades of the 20th century had worked on the physics of 
gas discharges, among them Rutherford, John S. Townsend, Owen W. Richardson, Karl 
Compton, Irving Langmuir and Peter Debye. In particular Townsend’s (Townsend 1901, 
1915) and Richardson’s (Richardson 1908, 1928) investigations of the different types of 
discharges were crucial for the further development of gas discharge physics. Clearly 
the peak of gas discharge physics coincided with the development of atomic physics and 
quantum mechanics, i.e. with the quantum theory of the atom. 

The first three decades of the 20th century saw the physics of gas discharge phenomena 
blossoming. With the investigation of the various processes of the ionisation of gases, by 
either radiation or collisions, the understanding of atomic spectra and of the excitation 


PAPER_SSRv.tex; 10/04/2008; 19:52; p.2 


Atmospheric discharge physics 3 


processes achieved quite a high level of sophistication and maturity. At the same time the 
physics of gaseous discharges bifurcated into two quite separate disciplines, proper gas 
discharges and plasma physics. The notion of a plasma was introduced by Irving Lang- 
muir in 1927 (Mott-Smith 1971) reserving it to highly ionised gases, with only a small 
residual admixture of a non-ionised neutral component or even lacking such a component. 
Due to the dominance of electromagnetic interactions in a plasma, the physics of plasmas 
turned out to be completely different from that of neutral gases. Peter Debye and Erich 
Hiickel (Debye and Hiickel 1923) showed that electrical charges in a plasma are screened 
beyond a distance Ap = vu. / Wye, termed the Debye radius, where ve = Vthe 18 the electron 
thermal velocity, and wy. is the electron plasma frequency, the oscillation frequency of 
electrons around their equilibrium position (discovered by Langmuir 1923). 

The notion that the upper atmosphere of the Earth must be electrically conducting 
goes back to Carl Friedrich Gauss (1776-1855) who, expanding the Earth’s magnetic field 
in spherical harmonics, discovered that it contained a substantial component due to outer 
sources which he suspected to be currents somewhere at large altitude above the surface. It 
was only in 1902 that Oliver Heaviside (Nahin 1990, Appleton 1947, Watson-Watt 1950) 
proposed the existence of a high-altitude conducting layer that reflects radio waves. This 
was confirmed experimentally in 1924-1927 by Edward V. Appleton (Appleton 1947). 
The formation and nature of this layer was understood when Sydney Chapman (Chapman 
1931) realised that solar UV radiation was capable of ionising the upper atmospheric lay- 
ers at around 100-130 km altitude, being absorbed there. H. Kallmann (Kallmann 1953) 
ultimately proved that taking into account the full solar spectrum the ionisation effectively 
reached down to below 80 km altitude. With the advent of the space age it became clear 
that atmospheric constituents below 80 km can also be ionised by precipitating energetic 
particles from the Sun and the magnetosphere and by cosmic rays. These processes are 
reviewed elsewhere in this issue (Chapter 2 and Chaper 6). 

In the early 1930s several extended reviews were written on this subject, among them 
the famous accounts of Compton and Langmuir (Compton and Langmuir 1930, 1931). Af- 
ter the golden twenties and thirties the physics of proper gas discharges lost its importance 
in fundamental physics, becoming a branch of applied physics, while plasma physics — 
its fully ionized twin — shifted into centre stage. This was mainly for two reasons, the 
interest in fusion research and the advent of the space age. Space plasma physics became 
attractive to both the public and politicians. Still, interest in ionospheric physics and in 
the mechanism of lightning enabled the physics of gas discharges to stay alive throughout 
this period of relative drought. 

In the last twenty years gas discharges became more interesting again, in relation to 
natural transient luminous events (TLEs) (Franz et al. 1990, Sentman et al. 1995, Neubert 
2003) that had been discovered in the Earth’s upper atmosphere above thunderstorms. 
Such spectacular phenomena are sprites, electric jets, and elves (Fullekrug et al. 2006b); 
evidence for lightning and related electrical effects on other planets, in planetary rings 
and even in astrophysical objects, accumulated. In addition, the relation of TLEs to space 
weather phenomena has provided another impetus to their investigation. [The processes 
of the generation and separation of charges which are required before a discharge can be 
generated are reviewed in this issue (Roussel-Dupré et al. 2008, Yair 2008) and also form 
the content of Chapters 2 and 6. Some of the theory and the observation of such effects 
as blue jets and sprites are summarized in Chapter 7 (Mishin and Milikh 2008, Mika 
and Haldoupis 2008).] Because of the importance of such effects the physics of gaseous 
discharges is far from being exhausted, in particular in view of the newly discovered fine 
structuring of streamers (Ebert et al. 2006). 
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There is one particular aspect of this that in the future will certainly regain more 
attention. This is the problem of the energetic coupling between near-Earth space and 
the atmosphere. More generally, the question of how a planetary atmosphere couples to its 
space environment poses a problem that is of fundamental importance for climatology and 
the understanding of how energy is fed into processes that are important for the planetary 
climate from the outside, the interplanetary and magnetospheric plasmas through coupling 
via the ionosphere to the atmosphere. Sprites, elves and jets suggest that this is indeed an 
important problem which has only just been realised and has not yet been investigated 
to the extent that it deserves. There are extended recent reviews of these phenomena 
available (Fullekrug et al. 2006b, Pasko 2007, for instance). 


2. Parameters and Theory 


The electric discharge is the final step of a whole sequence of processes going on in a 
gas that is subject to ionisation and internal gas dynamics. Most of these processes have 
been or will in subsequent sections in this issue be described for Earth and other planets 
surrounded by an atmosphere or some equivalent to that. At the end of all these processes 
there is an electric field of strength E that is maintained in some way over a sufficiently 
long time, and positive and negative electrical charge layers that are separated from each 
other by a distance L. This electric field can be maintained only up to a certain strength 
|E| < |E|crit which, when exceeded, causes a spectacular breakdown of the field. The 
processes of of electric field build-up are related to the various processes of ionisation and 
fake up considerable space in this issue. Here we restrict ourselves to a very basic and by 
no means complete discussion of the breakdown processes. Several aspects of these will 
be given in great detail in Chapter 7. 

The atmosphere is considered to be collisional, with the collision frequency 1. being 
determined mainly by collisions between electrons and the various neutral components of 
the atmosphere. A more sophisticated theory must take into account the composition of 
the atmosphere including aerosols, various excitation and ionisation processes such as, for 
instance, cosmic ray ionisation or frictional ionisation, the presence of a magnetic field, 
the horizontal and vertical dynamics of the atmosphere, its water content, temperature 
distribution, and its altitudinal inhomogeneity. All these factors (including those we are 
not considering here) enter into the real discharge process, and some of them are consid- 
ered in great detail elsewhere in this issue (Roussel-Dupré et al. 2008). Here we merely 
mention them as hard-to-treat complications of a realistic atmospheric discharge process 
encountered, for instance, in a lightning discharge or the recently discovered high-altitude 
discharges in sprites, blue jets and giant jets. 

In the simplest model of a resistive atmosphere the collisions are treated classically 
(and we will not deviate from this in this introductory paper). This implies that 0, % 4.5 x 
10-18 m~? (Aguilar-Benitez et al. 1990, Bernshtam et al. 2000) is the constant classical 
two-body collisional cross section. In addition the neutral gas density obeys the baromet- 
ric law n,(h) = no exp[—mng(h)h/ kpT,(h)], where the subscript n means neutrals, h 
is the height above sea level, no the density at ground level, g(h) the gravitational acceler- 
ation at altitude h, m the mass, T(h) the temperature at h, and kg Boltzmann’s constant. 
By contrast in a fully ionised plasma the cross section is the Rutherford or (when averaged 
over the thermal distribution function) the Spitzer-Braginskii cross section (Baumjohann 
and Treumann 1996, Chapter 4) which strongly depends on plasma density and particle 
speed, thus implying a completely different physics. Since in the uppermost layers of the 
atmosphere and in particular in the ionosphere the gas is weakly ionised with an admixture 
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Figure 1. Left: Average altitude profiles of temperature T'(h) and mass density p(h) in the Earth’s 
atmosphere consisting of 80% N2 and 20% Ox. Right: Rough models of collisional mean free path 
Amfp based on the data on the left for neutral-neutral binary collisions with collision frequency 
Ve = Vnn(h), and the electron-neutral collision frequency for collisions between electrons of 100 
eV energy and neutrals as functions of altitude above ground level. The light line shows the altitude 
dependence of the binary collisional mean free path Amfp(h). In the middle atmosphere the mean 
free path and collision frequencies vary about exponentially with altitude. The dashed vertical line 
shows the electron cyclotron frequency wee ~ 4.5 x 10° Hz which is nearly constant over this 
height range. The estimate of the electron collision frequency neglects excitation and Coulomb 
interactions which should be included above roughly 10 eV electron energy for a more realistic 
model. This would lead to a substantial modification of vy. = Ven. 


of electrons and ions, the relevant cross section will be a mixture of the classical and 
Spitzer cross sections. In the neutral atmosphere on the other hand several other effects 
have to be considered which are related to the excitation cross sections of the different 
gas molecules, the molecular nature of the gas and the composition of the atmosphere at 
different altitudes, as well as the different ionisation energies. Moreover, the eventual 
presence of an electric field and the background convective wind motions of the gas 
as well as the effects of the external radiation change the conditions existing. All such 
effects are ignored here but for realistic cases should be taken into account in their relative 
importance. 
With the above notations the collision frequency becomes 


ng(h)h T,,(0 2hkpTi 
v.(h) = vp exp | | 4/ at Vo = N00e4] mn, e (1) 


On Earth, the temperature and density of air at ground level are roughly Tp ~ 300 K and 
no © 2 x 107° m~®. The latter value depends on what is assumed for the composition of 
air. In the above formula the height dependence of the gravitational acceleration g(h) ~ 
g(0)(1 — 2h/Rg) can safely be ignored since in the altitude range from 0 to 100 km 
the variation of the ratio of height to Earth radius Rg is < 1%. The temperature (see 
Figure 1) in this altitude range changes by a factor of T,,(h)/To < 3. Thus its change is 
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important only in the exponential factor. At zero level, the Earth’s surface, we have for 
the neutral-neutral collision frequency vp ~ 10° Hz. 

Figure | summarises the vertical profile of the Earth’s atmospheric temperature, mass 
density and binary collision frequency in a simplified way with the neutral-neutral colli- 
sion frequency Vv. = Vp» and electron-neutral collision frequency v. = Ve, calculated for 
an 80% nitrogen and 20% oxygen mass-dominated atmosphere. The reason for the latter 
being large is the proportionality of the collision frequency to the thermal velocity uv. of 
the electrons which are generated by ionisation and therefore have energies > a few eV. 
At higher electron energies as those they gain during a discharge this collision frequency 
becomes obsolete because other processes (excitation, Coulomb interaction) set on. Since 
the Earth is magnetized the collision frequency is, in principle, a tensor (Rycroft et al. 
2008, e.g.) giving rise to a tensorial electrical conductivity 0 = (o\ ,Op,0H), where op 
is the conductivity perpendicular to b = B/B, the unit vector parallel to Earth’s magnetic 
field B, and oy is the Hall conductivity perpendicular to both, the electric and magnetic 
fields. The expressions for the components of o are o = €? Ne /MeVe = E0W ee /Ve 
parallel to b, and 


2 


V Ve W, 
Cc (63 ce 
oP =z || OH>=—- a 3 4l (2) 
2 2 7 2 2 : 
Ve st Wee UY zn Wee 


where wee = eB/me and Wye = e\/n/eome are the electron cyclotron and plasma 
frequencies, respectively, and € is the dielectric constant of free space. Below the meso- 
sphere these expressions simplify to cp = O|,0H ~~ —Wee /Ve < 1. Hence the at- 
mosphere is practically isotropic with zero Hall conductivity. As Figure | suggests, this 
changes at altitudes above 90 km at the top of the mesosphere. Here the collision fre- 
quency drops below the electron cyclotron frequency (the vertical line in the figure) 
which has a weak altitude dependence only according to the weak variation of the Earth’s 
magnetic field B ~ Bo(1—3h/Re») with height h through the atmosphere; Bo is the value 
at the Earth’s surface on a particular geomagnetic field line. The full conductivity tensor 
comes into play, and this has consequences for the development of electric discharges at 
these altitudes. The conductivity oj is low in the atmosphere meaning that the atmosphere 
is a poor conductor. Under discharge conditions this may change, though. In contrast, 
the ions are practically collisionally isotropic because vy, >> w,; by far exceeds the ion 
cyclotron frequency w.; = Z(me/m;) Wee up to the altitudes of the ionosphere (Z is the 
ionic charge number, and m, the ion mass). 

Introducing the more detailed papers collected in Chapter 7 on the most interesting 
special cases of atmospheric discharges in the upper atmosphere and their effects on radio 
wave propagation, we will, in the following, briefly review the types of electrical dis- 
charges that may be relevant under gaseous atmospheric conditions. We explicitly exclude 
here other important discharges such as those in charged dust and sand storms which are 
considered in detail elsewhere in this volume (Renno and Kok 2008). On Earth they occur 
in volcanic eruptions which inject large amounts of dust and ashes into the atmosphere 
where they trigger lightning. They are also known to occasionally occur in sand storms in 
the big terrestrial sand deserts. In the atmospheres of the dusty planets Mars and Venus 
they become important if not dominant (Melnik and Parrot 1998, Farrell et al. 1999). 


3. Different Types of Discharges 


In order to ignite an electric discharge in an otherwise neutral gas like the atmosphere, 
electrons have to be set free by some process. Subsequently an electric field must be 
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Figure 2. Two types of electric discharges. Left: A typical arc discharge, and Right: a glow discharge 
in an evacuated glass tube. Glow discharges are quiet and low temperature which is in contrast to the 
noisy and unsteady spark discharges. Because of this reason they can be used for illumination. Arc 
discharges, on the other hand, develop high temperatures and find technical application in welding. 


present that can accelerate electrons to energies far above the ionisation energy fast enough 
in order to cause an avalanche of newly generated electrons by collisional (or other kinds 
of) ionisation. For this to happen the electric field must be strong, which implies that by 
some external force the initially present negatively and positively charged layers must 
become separated over a sufficiently large distance L. Processes capable of providing 
such initial charge distributions are discussed elsewhere in this issue. One of the basic 
processes is the weak ionisation of the atmosphere produced by the continuous inflow 
of highly energetic cosmic rays into the atmosphere and the cascades of nuclear fission 
products and elementary particles which are produced by them. Charge layer separation 
in a thunderstorm is a complicated process closely related to the physics of clouds and 
their internal dynamics. The reader is referred to the respective sections of this volume 
for information. 


3.1. TOWNSEND, GLOW, AND SPARK DISCHARGES 


Once sufficient charge has accumulated to form a charge layer and oppositely charged 
layers have become separated from each other, a macroscopic electric field is built up. 
As long as this field is weaker than some threshold field E.rit, discharges go on only 
slowly due to recombination and may be balanced by newly created ionisation. However, 
when the electric field exceeds the threshold, a violent discharge process sets in, causing 
breakdown and the electric field is shorted out. In laboratory gas discharge physics one 
distinguishes a number of discharges. Two of them, which are accompanied by light 
emission, are shown in Figure 2. To the left is the typical blue arc discharge between 
cathode and anode, and to the right is a typical red glow discharge observed in a glass 
tube at low pressure. 

In addition there are dark (or Townsend) discharges which, as the name implies, are 
not accompanied by visible light emission. Moreover, one distinguishes non-stationary 
discharges in connection with spark discharges. In Townsend dark discharges the current 
flow is rather weak and the current densities are low; the voltage remains constant over the 
time of the discharge. When the current strength increase the Townsend discharge makes 
the transition to a glow discharge. Now the energy of the flowing electrons is large enough 
to excite an atom collisionally, which emits light from the infrared to the visible and up 
to the ultraviolet. Arc discharges acquire large current densities, strong electric fields, 
and may emit radiation deep into the ultraviolet and even weak X-rays. These processes 
have been described in textbooks (Lieberman and Lichtenberg 1994, Raizer 1997, for 
instance). In atmospheric physics of most interest are discharges of the kind which lead 
to breakdown as, for instance, when lightning discharges are generated. 
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Figure 3. Left: Avalanche generation in electron collisions with neutrals in presence of a suffi- 
ciently strong electric field (the Townsend process). Right: The resulting charge distribution in the 
avalanche. The head of the avalanche consists of fast and dense electrons. An extended wake of slow 
ions is left behind which try to follow the electrons but are accelerated in the opposite direction by 
the externally applied electric field which, however, is partially cancelled inside the ion cloud wake. 


3.2. ELECTRON AVALANCHES 


For the purposes of atmospheric discharges the formation of current carrying electron 
avalanches is most important, as shown in Figure 3. The continuous slow discharges 
of electric fields produced by the persistent inflow of cosmic rays into the atmosphere 
provide the seed population of electrons for the ignition of an avalanching discharge. In 
an avalanche the electron number density increases exponentially since dn-/dt = ajne, 
the production of electrons being proportional to the existing electron number density and 
with the proportionality factor being the Townsend ionisation rate coefficient a; (Raizer 
1997). This relation has the solution n-(t) = ne(0) exp[f, dt'a;(t’)]. The exponential 
factor indicates the temporal multiplication rate of the avalanche, which also corresponds 
to the spatial multiplication rate over the distance of the evolution of the avalanche. 

Avalanches have a typical structure with a negative head consisting of the fast electrons 
and an extended positively charged tail containing the slow freshly produced ions which 
tend to retard the advancement of electrons but are kept back by the external electric 
field that accelerates them in the opposite direction. It is clear that in such a process the 
electric field must become strong enough in order to overcome recombination collisions 
and accelerate electrons in one mean free path length (see Figure 1) up to ionisation 
energies. The latter are of the order of a eV, for hydrogen typically Ein. ~ 13.6 eV. 
The accelerated electrons must be able to both run away and ionise. From this reasoning 
follows a simple condition for the necessary strength of the external field E for igniting a 
gaseous electric discharge 


|E| z Eerit = 2(Eion/€Amfp)- (3) 


At 10 km altitude this field is quite large, amounting to > (2 — 3) x 10° Vm7}, but 
decreases with altitude deep into the mesosphere (~ 70 km altitude) by roughly four 
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orders of magnitude (see Figure 1) to become only a few 100 V m7! there. Hence, in 
the higher atmosphere considerably smaller electric field strengths are required in order 
to ignite an electric discharge than in the troposphere. The problem therefore consists less 
in exceeding the threshold at these altitudes than in generating charged layers with an 
electric field between them which requires vertical transport. But once this is achieved 
at lower altitudes breakdown at higher altitudes becomes possible. In this respect it is 
noteworthy that sprites and jets have been observed to be closely related to low altitude 
thunderstorm discharges starting from the tops of thunderstorm clouds at < 10 km. 

The high threshold electric field must in fact be exceeded at low altitudes of a few 
km, which poses a serious problem. A solution has been proposed when it was theoret- 
ically realised and proved by kinetic simulations (Gurevich et al. 1992, Roussel-Dupré 
et al. 1994, Lehtinen et al. 1997, Babich et al. 1998) that a proper kinetic description of 
the avalanching process should take into account not only the generation of secondary 
electrons but also the production of high energy run-away electrons. This causes what 
these authors call “runaway breakdown” (RBD). A complete review of the theory can be 
found elsewhere (Gurevich and Zybin 2001). The idea is that a small number of fast seed 
electrons initially reach high energy (= 10 keV) in the electric field. At these energies 
the electrons are fast enough to ignore the neutral state of matter. They take the outer 
shell valence electrons and the nuclei of the air molecules as free particles. For them the 
ordinary collisional cross section @¢,, for ionisation is replaced with the energy dependent 
Rutherford-Coulomb cross section og when passing through the interior of the neutral 
molecules. This drastically increases the efficiency of ionisation since the Coulomb cross 
section is inversely proportional to the square of the particle energy; it decreases with 
energy which allows for energetic particles after collision to “run away” in the field, and 
obey a power law distribution function f(€) « €~1! at energies below some threshold 
energy above which they lose their motional energy by radiation. It has been shown 
(Gurevich et al. 1992) that this self-consistent power law distribution is a function of 
the ratio of the electric to breakdown-threshold electric fields. 

The RBD threshold electric field as a function of altitude h is proportional to the neutral 
atmospheric density 


Mn (h) Mngh 5 V 
E| = Eri ~ 0.2 | ——| = 2.0 € ——_ 10° —. 4 
|E| Zz t,RBD | am exp | kT, (h) x ma (4) 
This value is one order of magnitude less than the value obtained from Eq. (3). Due to the 
dependence of n,,(/) on altitude, it decreases exponentially with h, and at higher altitudes 
becomes very low indeed. The Townsend coefficient now becomes a function of the ratio 
of electric field to critical electric field (Gurevich et al. 1992) 


O;RBD & 0.07(E/Eorit.RBD)?; (5) 


showing the pronounced nonlinearity of the self-feeding process, i.e. avalanche formation. 
The numerical factor 0.2 in Eq. (4) has been determined from numerical simulations of 
runaway discharges at different pressures. 


4. Secondary Effects 


The avalanches produced in these spark discharges in the atmosphere cause a number of 
secondary effects. Sparks are bright, emitting light due to the excitation of both atoms and 
ions which emit in various lines of the electromagnetic spectrum. Moreover, space charges 
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that are built up over macroscopic vertical and horizontal length scales by the dynamics 
of the neutral atmosphere due to moving the charges act back on the dynamics, braking 
it due to the large electrostatic forces that develop when the charges become sufficiently 
large. Finally, the ionisation of the atmosphere affects the atmospheric chemistry. Other 
effects are the excitation of several kinds of waves, the generation of radio waves and 
modification of propagation conditions. These have been observed to occur during and 
in the wake of thunderstorms and in coincidence with sprites and jets (Bosinger and 
Shalimov 2008, Mika and Haldoupis 2008, Simoes et al. 2008). Low frequency elec- 
tromagnetic radiation (VLF, ELF, ULF, Schumann resonances) which accompanies TLEs 
is well understood (Rodger 1999, Reising et al. 1996, Cummer and Inan 1997, Fiillekrug 
et al. 2001, Fiillekrug et al. 1998, Fiillekrug and Reising 1998, Fiillekrug et al. 2006a). 
In the following we list some secondary microscopic processes in avalanches but are 
less familiar. Some of them might affect the physics of the discharge and might result 
in observable effects. 


4.1. WAVES 


Ion sound. The simplest wave type is the mechanical distortion of the atmosphere 
during discharges which cause thunder and mechanical oscillations of the atmosphere. 
During such disturbances the ionic component of the atmosphere undergoes similar distor- 
tions which may result in low frequency (0 < wis /Wpe<1/Me/mM;) Sonic waves Wis = keis 
of speed vj; < Cis ~ \/Me/M; Ve ~ 40km s—+, where ve, v; are the thermal velocities of 
electrons and ions, respectively. 

It is usually believed that these waves are strongly damped due to the high collision 
frequency at low altitudes (see Figure 1). Landau damping can be neglected because of 
the low ion temperature. At higher altitudes than at the top of the mesopause ion sound 
turbulence should accompany any discharge. We must compare the mean free path of 
electrons with the Debye length Ap=ve/wpe, where the plasma frequency wpe ~ 60,/Ne 
s~+, and the electron density is measured in m~°. The electron density produced in a 
discharge is some small altitude-dependent fraction ¢(h) <1 of the neutral density n,,(h). 
The ratio of the Debye length to the mean free path Amp = 1/nn(h)ve is AD/Amftp ~ 
10-?,/(kpTe/e)(mn(h)/noG(h)) < 1 for any finite electron temperature of a few 10 eV, 
and ¢(h) < m»(h)/no. Hence, the plasma properties dominate inside the avalanche. 

Moreover, the electrons being of energy 210 eV, are much hotter than the ions, T,. >> 
T;. They carry the current 7 =—en,.Va, where Vq is the current drift velocity. This makes 
excitation of ion sound waves probable if only the growth rate y;, > v- exceeds the col- 
lisional electron damping rate. In an avalanche the latter is given by the Spitzer collision 
frequency vc. This condition, which can be written 1 > Vg/ve > 200/Np, is always 
satisfied, then. With Np = (no: >> 1 the number of electrons in a Debye sphere, its 
right-hand side is a small number. The ion-sound waves cause the avalanche to develop 
magnetic field-aligned striations (Davidson 1972, Treumann and Baumjohann 1997, e.g., 
for the discussion of the instabilities and their nonlinear evolution). 


Buneman two-stream instability. For Va > ve the avalanche plasma is Buneman two- 
stream unstable (Buneman 1958). RBD discharges with their comparably high electron 
density and electron kinetic energy (Gurevich and Zybin 2001) are candidates for the 
two-stream instability. The two-stream instability is riding on the beam, i.e. it propagates 
at about the same velocity as the beam, being for long time in resonance with the current 
carrying electrons and thus giving rise to high wave intensities. This instability saturates 
by formation of localised electron (and also ion) holes, which correspond to locally (a 
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Figure 4. Electron phase space plots at three successive simulation times (measured in plasma times 
Wye) of the formation of electron holes in a two-stream unstable plasma. Instability is forced by 
the presence of a density dip (potential wall for electrons) in the centre of the simulation box which 
interacts with an injected electron beam current (after Newman et al. 2002). The box has one space 
(measured in Debye length) and one velocity coordinate (measured in electron thermal speeds ve. 
Hole formation is seen on the right. The holes move to the right with beam velocity. They trap 
electrons and cause strong electron heating of the trapped electrons. In addition acceleration of a 
narrow (cold) beam is observed which escapes from the hole region at about 7 times the initial 
velocity. 


few 10 Debye lengths long) very strong electric fields. This has been demonstrated in 
numerical simulations (Newman et al. 2002). [The existence of holes has been confirmed 
in the auroral plasma (Carlson et al. 1998, Ergun et al. 1998, Elphic et al. 1998).] De- 
pending on its polarity, the localised electric field traps electrons or ions. The trapped 
avalanche electrons bounce back and forth in the holes and become violently heated, 
dissipating a substantial part of the beam current and electric field energy. Conversely, 
passing (untrapped) electrons are accelerated further and are at the same time cooled to 
low temperatures, as seen in Figure 4. 


Phase space hole effects. The holes resulting from the two-stream instability are re- 
gions of high electric field pressure which affects the equilibrium between the avalanche 
and its neutral environment. Force balance requires V,P, = —VePe—- VePer, where 
P, =1nkpT,, is the atmospheric pressure, P. = n-kpT, is the electron pressure in the 
avalanche, and Pg = }€9|E|?{0|[we(w)] /Ow} is the field pressure. The electron and field 
pressures together in the avalanche act to expel the neutral gas from the avalanche region. 
But the forces on the right and left act on different scales. The scale height L,, of the 
pressure is several kilometers, while the scales L. ~ Lg of the holes are only a few 
10 Debye lengths long. The local force balance thus becomes 


Ly 


P, = 
10\p 


(P. + Pr). 
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Moreover, P,, ~3 x 104n»(h)/no Jm~3. Assume that the electrons have an energy of 10 
keV (just a fraction of the energy assumed in Gurevich’s RBD mechanism). Then, P. ~ 
10~!°n, Jm~®. To calculate the electric field pressure one needs the wave dielectric con- 
stant €(w). For the two-stream instability we have O[we(w)|/Ow~w%,/(0.03 wpe)? ~ 10°. 
Thus, Pg ~ 5x 10~°|E|? Jm~%. The Debye length is \p ~ (0.7/\/ne) m. The balance 
equation becomes 


1041p (h) ~ 102L,,(km)n2 (107 n, + 5 x 10-9] E[?). 
This condition can be satisfied with moderate avalanche densities n. < 10/9? m~? and 
electric fields |E| ~a few Vm7~!. Stronger fields and higher avalanche densities will 
completely blow the neutral atmosphere away from the region of the hole. This crude 
estimate shows that the localisation of the electric field in electron holes and trapping of 
electrons as a consequence of the two-stream instability causes a substantial local dilution 
of the neutral atmospheric gas in the holes in the avalanche. It causes a filamentation of 
the avalanche into striations and narrow packets of trapped electrons and strong fields. 
Fine structuring has been observed in streamers (Ebert et al. 2006) though this impor- 
tant observation is interpreted differently. Ebert et al. (2006) describe their observations 
as a multiscale structure of streamers occurring mainly in the head of the streamers. 
We should, however, note that the holes ride on the beam and move with about beam 
velocity. Thus the effect described here also happens mainly in the narrow head layer of 
the avalanche where the electron density, energy and velocity are highest. 


Local whistlers. Avalanches can serve as sources of locally excited whistlers. This is 
obvious from considering the trapped electron distribution in Figure 4. The electrons that 
are trapped in the electron hole potentials bounce back and forth along the magnetic field. 
Due to their gyrational motion in the geomagnetic field they possess magnetic moments 
[te =T.1 /B (Te. is the electron temperature perpendicular to the geomagnetic field). As 
a consequence of the bounce motion of the trapped electrons, the electron holes become 
unstable (Ergun et al. 1998) with respect to whistlers of frequency w.;§<w<<w,-. While 
moving with the parallel hole velocity along the magnetic field, they radiate whistlers act- 
ing like point sources. In a spatially varying magnetic field this parallel motion gives rise 
to a saucer like frequency-time appearance of the whistlers. [Such whistler-saucer emis- 
sions were indeed observed first under the conditions in the auroral ionosphere (Ergun 
et al. 1998). Similar phenomena during thunderstorms have been reported (Parrot et al. 
2008) from the low-altitude spacecraft DEMETER.] Their propagation direction is oblique 
to the geomagnetic field at the whistler-resonance cone, with frequency near the lower- 
hybrid frequency wp, ~ wp; /(1+w7./w2.) 2, In the strong geomagnetic field this is close 
to the ion plasma frequency w,;. Their frequency-time dependence maps the avalanche 
plasma density along the whistler path. Propagating obliquely, these short-perpendicular 
long-parallel wavelength whistlers are trapped in the avalanche. Being multiply reflected 
from its boundaries, the whistlers ultimately escape on the top-side with Alfvén veloc- 
ity V4 = B/\/Hommne mostly in perpendicular direction to the magnetic field thereby 
spreading over a large hollow horizontal area of narrow opening angle, centred around 
the source before being absorbed. Excitation of the mesospheric or exospheric gas might 
lead to emissions over the whole area of these whistlers resembling the Elve emission. 
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Figure 5. Left: The L+L’ — T plasma wave process for generation of electromagnetic radiation 
T. The large parallel wave numbers kj of the counterstreaming Langmuir waves cancel, and a 
transverse wave of frequency w3 = w1 +W2 & 2 Wye (with very small perpendicular wave number 
k3~~k1 <hky and therefore large phase and group velocities ~ c) is emitted. Right: Scattering of 
an L-wave by an ion i as a process in which an L’-wave is generated. i* is the excited ion which is 
left over after the ‘collision’ with the L-wave. 
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4.2. RADIATION 


Several types of electromagnetic radiation have been observed in relation to atmospheric 
electric discharges, both natural and artificially triggered (Dwyer et al. 2003) lightning, in 
addition to the above mentioned spherics and the whistlers propagating in the ionosphere 
and magnetosphere caused by them. Lightning discharges have been observed to be ac- 
companied by radiation in the optical, UV, X rays and even up to the gamma ray energy 
range and, also in the ELF and ULF radio wave ranges (Simoes et al. 2008, Bosinger and 
Shalimov 2008, Mika and Haldoupis 2008). Moreover, spacecraft observations of other 
planets (Jupiter, Saturn) provide evidence of short duration radio emissions in connection 
with lightning on those planets (Fischer et al. 2008). 


Gamma- and X-rays (Bremsstrahlung). Referring to the above discussion of the 
general physics of breakdowns it can be expected that radiation is emitted from breakdown 
discharges in the optical to X-ray ranges. An early prediction (Wilson 1925, 1956) based 
on the Townsend mechanism that lightning should be accompanied by the generation of 
energetic electrons and possibly by the emission of X-rays has been confirmed by the 
more recent BATSE experiment aboard the Compton Gamma Ray Observatory (Fishman 
et al. 1994) and by the Alexis spacecraft observations (Blakeslee et al. 1989, Holden et 
al. 1995). Gamma radiation is believed to be caused by the energetic electron component 
generated in the runaway phase of the discharge when avalanches of electrons are formed. 
In the light of the above discussion on holes it is the fast cold high energy electron beam 
that is responsible for their emission. The radiation mechanism has been investigated 
theoretically (Gurevich et al. 1992, Roussel-Dupré et al. 1994). Since very high energy 
electrons are required to generate these Gamma ray bursts the observations completely 
falsify Wilson’s mechanism while they are in excellent agreement with the RBD mecha- 
nism, providing a strong argument in favour of runaway breakdown and lowering of the 
avalanching threshold. The gamma rays observed are not caused in a nuclear interaction, 
and no gamma lines have been detected yet. They form the high energy tail of the energetic 
X ray emission resulting from free-free (‘bremsstrahlung’) radiation of the most energetic 
runaway electrons. The absence of lines is an indication for the existence of an upper 
threshold of the accelerating electric field. 


Radio waves: Plasma processes. Any emitted high-frequency radio radiation can 
hardly be via the synchrotron process (Rybicki and Lightman 1979) as this requires much 
higher emission measures, i.e. large volumes occupied by energetic electrons in addition 
to much stronger magnetic fields, than those which are available. The only imaginable 
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mechanisms are plasma mechanisms like the head-on interaction of beam generated Lang- 
muir (L) waves due to a process L+L’ — T (where T is the transverse long-wavelength 
radio wave). This process is shown schematically in Figure 5. 

Generated by the cold fast electron beam, L-waves propagate in the direction of the 
beam. The backward moving L’-wave can be generated by backscattering of the L-wave 
by the slow cold ion component according to the process L+i—L’+i*. Here i is the involved 
ion, L’ the backscattered Langmuir wave, the wave vector direction of which is reversed, 
and i* is the excited ion (which emits light). L-waves can also interact with low frequency 
ion sound waves (IS) according to the reaction L+IS —T, following the same scheme as 
shown in Figure 5. This process has a lower efficiency than the L+L’ process, though. In the 
former case the frequency of the emitted radio wave is close to twice the plasma frequency 
w3(L+L) ~ 2Wye. Because the frequency wis << Wye of the S-wave is much smaller than 
the electron plasma frequency, the process L+IS emits radiation of frequency near the 
plasma frequency: w3(L+IS) ~ w»-. For avalanche beams of density n. ~ 10'°m~? this 
frequency is about ~1 MHz. 

Excitation of electron-acoustic waves (EA) of frequency close to a fraction of wpe 
can also contribute to HF-radio emission. The reaction equation is L+EA—T, and the 
frequency is above but close to Wye, typically w3(L+EA)S 1.5 wpe. All these mechanisms 
depend crucially on the presence of a hot plasma electron component, of cool ions, and 
only a small admixture of cold background electrons. These requirements imply that a 
substantial fraction of the neutral atmospheric components has been locally expelled from 
the avalanche region. 


Radio waves: Electron-cyclotron maser. Very intense HF-radio waves can be gen- 
erated by the electron-cyclotron maser (Treumann 2006, for a recent review). This is 
an extraordinarily efficient process. Its requirements are a strong field-aligned electric 
potential difference, parallel current flow and a strong geomagnetic field in addition to 
a negligible component of cold electrons and neutrals. It also requires some pitch angle 
scattering of the electron distribution such that the electrons have a substantial velocity 
component perpendicular to the geomagnetic field. Again, this mechanism can be realised 
only if the neutral atmospheric gas is expelled from the avalanching electron cloud. The 
electron-cyclotron maser is a direct process in which no other catalysing waves are in- 
volved. It generates long wavelength (Accem~ few 100 m) electromagnetic radiation at the 
local electron cyclotron frequency Weem ~ Wee, Which in the atmosphere is in the MHz 
range. The bandwidth of the emitted radiation depends on the electron beam temperature 
ranging from very narrow (< W,-) up to a substantial part of the electron cyclotron 
frequency. 
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